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ABSTRACT 

We employ 76 type la supernovae with optical spectrophotometry within 2.5 days of B-band maximum light obtained by the Nearby 
Supernova Factory to derive the impact of Si and Ca features on supernovae intrinsic luminosity and determine a dust reddening law. 
We use the equivalent width of Si n /14131 in place of light curve stretch to account for first-order intrinsic luminosity variability. The 
resultant empirical spectral reddening law exhibits strong features associated with Ca n and Si n /16355. After applying a correction 
based on the Ca n H&K equivalent width we find a reddening law consistent with a Cardelli extinction law. Using the same input data, 
we compare this result to synthetic rest-frame UBVRI-like photometry in order to mimic literature observations. After corrections for 
signatures correlated with Si n /14131 and Ca n H&K equivalent widths, and introducing an empirical correlation between colors, we 
determine the dust component in each band. We find a value of the total-to-selective extinction ratio, R v = 2.8 ± 0.3. This agrees with 
the Milky Way value, in contrast to the low Ry values found in most previous analyses. This result suggests that the long-standing 
controversy in interpreting SN la colors and their compatibility with a classical extinction law, critical to their use as cosmological 
probes, can be explained by the treatment of the dispersion in colors, and by the variability of features apparent in SN la spectra. 
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1. Introduction 

Type la supernovae (SNe la) luminosity distances are measured 
via the standardization of their light curves using brightness- 
width (stretch, X\, Ami 5 ) and color corrections (Phillips 1993; 
Tripp 1998; Guy et al. 2007; Jha et al. 2007). While the deter- 
mination of the intrinsic dispersion related to light curve shape 
is subject to small differences between fitters, the manner in 
which color is linked to dust is still controversial, as it may 
be affected by additional yet unidentified intrinsic dispersion. 
Whereas earlier work used the total-to-selective extinction ra- 
tio of the Milky Way, Ry-3.1, direct estimates from supernovae 
(SNe) Hubble diagram fits lead to lower values, from Ry-l.l 
to R v =2.5 (Hicken et al. 2009; Tripp 1998; Wang et al. 2009). 
While the derivation of this value is subject to assumptions about 
the natural color dispersion of SNe (Freedman et al. 2009; Guy 



et al. 2010), the reason for a difference between SNe data and 
the Milky Way average result had remained unknown. 

Equivalent widths are good spectral indicators for addressing 
this issue as they probe the intrinsic variability of SNe la and 
by construction they have little dependence on extinction due 
to their narrow wavelength baseline. Arsenijevic et al. (2008) 
showed a strong correlation between the equivalent width of the 
Siii/14131 feature and the SALT2 x\ width parameter (Guy etal. 
2007), and Bronder et al. (2008) showed its correlation with M B . 
Walker et al. (2010) used it to standardize the Hubble diagram, 
but were hindered from drawing firm conclusions by the quality 
of the low redshift data. 

In this work, we take advantage of the Nearby Supernova 
Factory (SNfactory) spectrophotometric sample to revisit these 
conclusions, using both spectral data and derived UBVRI-like 
synthetic photometry. We present in Section 2 the SNe la sam- 
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pie and the definition of the Si n A4 1 3 1 and Ca n H&K equivalent 
widths, which will be used in Section 3.1 to correct the Hubble 
residuals. These corrected magnitudes are used to derive the rel- 
ative absorption in each wavelength band, 6A A . The correlations 
between the 5A A from SN to SN across different bands provide 
the reddening law, as described in Section 3.2, as well as the 
dispersion matrix between bands. In Section 4 it is shown that 
the resulting reddening law agrees with a Cardelli extinction law 
(CCM, Cardelli et al. 1989; O'Donnell 1994). Our reddening 
law has a value of Ry that agrees with the Milky-Way value of 
3.1, when the proper dispersion matrix is used. We then discuss 
these results in Section 5 and conclude in Section 6. 



2. Data-set and derived quantities 

This analysis uses flux calibrated spectra of 76 SNe la ob- 
tained by the SNfactory collaboration with its SNIFS instrument 
(Aldering et al. 2002) on the University of Hawaii 2.2-meter tele- 
scope on Mauna Kea. This subset is selected in the same way as 
in Bailey et al. (2009), using only SNe having a measured spec- 
trum within 2.5 days of B-band maximum, but with an enlarged 
data-set with a redshift range of 0.007 < z < 0.09. The SALT2 
X\ and c parameters (Fig.l), as well as the spectra phases with 
respect to the B-band maximum light, are derived using fits of 
the full light curves in three observer frame top hat bandpasses 
corresponding approximately to BVR. 

After flux calibration, host galaxy subtraction and correction 
for Milky Way extinction, the 76 spectra are transformed into 
rest frame. For the spectral analysis, they are rebinned with a 
resolution of 1500 km s _1 . In addition, synthetic magnitudes are 
derived from the spectra after photon integration in a set of five 
UBVRI-like top-hat bandpasses with a constant resolution over 
the whole spectral range (3276 - 8635 A, see Fig. 3 and caption). 
The uncorrected Hubble residuals 5M(A) - M(A) - (M(A)} are 
independently computed for each band of mean wavelength A, 
relative to a flat ACDM universe with Qm = 0.28 and Hq = 
70 km s Mpc~' . For the spectral analysis, the Hubble residuals 
are corrected for phase dependence using linear interpolation. 
For broad bands, we instead interpolate the magnitudes to the 
date of B maximum using the light curve shape of each SN la as 
defined by the SALT2 x\ parameter, but with each band's peak 
magnitude fitted separately. The data are presented in the online 
Table 1 . Error on Hubble residuals include statistical, calibration 
(~ 2% for peak luminosity) and redshift uncertainties. They are 
correlated between bands with a correlation coefficient varying 
from 0.64 to 0.98. 

The equivalent widths ew Si and EW Ca corresponding to the 
Si ii /14131 and the Ca n H&K features are computed in the same 
way as in Bronder et al. (2008, eq 1 ). The errors are derived using 
a Monte Carlo procedure taking into account photon noise and 
the impact of the method used to select the feature boundaries. 
ew Si and EW Ca are insensitive to extinction by dust, changing by 
less than 1% when adding an artificial reddening with a CCM 
law with R v = 3.1 and E(B - V) = 0.5. 



3. Method 

3.1. Intrinsic corrections and absorption measurements 

The Hubble residuals exhibit a dependency on observables such 
as ew Si and EW Ca which are uncorrected with dust extinction. 
As shown in Fig. 1 as an example, the 6M B dependence on 
ew Si exhibits a linear behaviour, with an asymmetrical magni- 
tude dispersion attributed to extinction and remaining intrinsic 




Fig. 1. Correlations of ew Si , SALT2 X\ (top left) and c (top right) pa- 
rameters and measured peak absolute magnitude up to a constant term, 
6M B (bottom). The open circles in the bottom panel are the data points 
excluded from the fit, which is displayed as a solid line, is equivalent 
for the B band to the red curve shown in Fig. 3a. 

variability. Similar dependence with equivalent widths is found 
for other bands. We may thus model the Hubble residuals for a 
given SN, /, as a sum of intrinsic and dust components, making 
various assumptions about the number of spectral energy distri- 
bution (SED) correction vectors, s A , from none (Eq. la) to two 
(Eq.lc): 



Ewfsf 

1 A 

EwfVf 
' A 



+ Ewps^ + SAff* 



(la) 
(lb) 
(lc) 



After finding and the <5A* by a least squares minimization 
over all SNe from Eq. lb, with an asymmetrical 2cr clipping, s^ 1 
is kept constant and s^ a and the <5A** Ca are then computed in the 

same manner. An example of the fit using ew Si for the B band 
is given is Fig. 1 (solid line), s^ being the correction applied to 
find the SA s j.. 



3.2. The empirical reddening law 

As we expect the relation between the 5A A to be linear for dust 
extinction, as shown in Fig. 2 for the U and V bands, we model 
the empirical reddening law as: 

6A AJ = y A 6A* VJ + 77, (2) 

where 5A Ai are the measured values, the slopes y A (- dA A /dA* v ) 
the reddening law coefficients, 5A* vj the fitted relative extinction 
for the SN i, and rj A a free zero point. 6A* vj , y A and r\ A are ob- 
tained by ax 2 fit using the full wavelength covariance matrix C, 
of the 6A A j measurements. However, the data are more dispersed 
than their measurement error (Fig. 2), and another source of dis- 
persion must be introduced. We adopt an iterative approach by 
using the fit residuals, r A j = 6A A j - (y A 6A* vj + rj A ), to determine 
the covariance remaining after accounting for the measurement 
error covariance, C AlAl j. This empirical covariance matrix, D, is 
given by: 



D 



A,A 2 



1 N 

~ N^j ( rXui rAl ' i ~ CaiA ^' 



i) 



(3) 
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Fig. 2. Relation between SAy and 6A V after ew Si correction (triangles 
up) and after ew Sj and EW Ca corrections (triangles down). The SAu are 
displayed with an added arbitrary constant. The ellipses represent the 
full measured covariance matrix between the two bands. They are en- 
larged by the additional subtraction error after EW Ca correction. The re- 
sult of the fits are also displayed. The y^ 1 and y^ 1+Ca respectively corre- 
spond to the ones in Fig. 3b and Fig. 3c for the U band. 
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Fig. 3. Black: Reddening law presented as y A = 8A A /dA* v as a function 
of wavelength. Filled circles correspond to the results obtained using 
the UBVRI-like bands, curves are for the spectral analysis. Red: Linear 
slope s A (mag/A) of equivalent widths versus SM A . Dotted lines: CCM 
law fit corresponding to the broad bands analysis. Panel a: 6M A cor- 
rected only for the phase dependence (eq. la). Panel b: SM A corrected 
for phase and ew Si (eq. lb). Panel c: 6M A corrected for phase, ew Si and 
EW Ca (eq. lc). The vertical dotted lines represent the UBVRI-like bands 
boundaries. The shaded vertical bands represent the Si and Ca domain. 
The shaded band around the curves are the statistical errors. 



where N is the number of SNe. In the next iteration, the total 
covariance matrix is given by C, + D. We have checked that the 
converged matrix does not depend on initial conditions. 



4. Results 

4.1. ew Si and Ew Ca impacts on the derived extinction law 

Results for the SED correction vector, s^, and the reddening law, 
jx, are presented in Fig. 3 for different assumptions about the 
number of intrinsic components. If SNe were perfect standard 
candles affected only by dust as assumed by Eq. la and Fig. 3a, 
the empirical reddening law y9 would be a CCM-like law with an 
average Ry for our galaxy sample. However, y° A clearly exhibits 
small-scale SN-like features. These features correlate strongly 
with some of the features in the ew Si correction spectrum, s s \ 
derived via Eq. lb and illustrated in red in Fig. 3a. 

The reddening law obtained after ew Si correction (Fig. 3b) is 
closer to a CCM law, except in the Ca n H&K and IR triplet, and 
the Si ii -16355 region. This indicates the presence of a second 
source of intrinsic variability. We select EW Ca to trace a second 
spectral correction vector, s Ca (Eq. lc), since Ca is clearly a ma- 

jor contributor to the observed variability and EW Ca and ew Si also 
happen to be uncorrelated (p = 0.06±0.12). As shown in Fig. 3b, 



s^" does a good job of reproducing the shape of the deviation of 



y\' relative to the CCM law. 

The mean reddening law, -y^' +Ca , obtained after the addi- 
tional correction by EW Ca (Fig. 3c) is a much smoother curve 
with small residual features, and agrees well with a CCM ex- 
tinction law. Thus it appears that these two components can ac- 
count for SN la spectral variations at optical wavelengths. Any 
intrinsic component that might remain would have to be largely 
uncorrelated with SN spectral features fixed in wavelength, as 
well as being coincidentally compatible with a CCM law. 

4.2. R v determination 

We apply the same treatment as above to our UBVRI-like syn- 
thetic photometric bands, and find agreement with the spec- 
tral analysis, as shown by the black points in the three pan- 
els in Fig. 3. After the ew Si correction (Eq. lb, Fig. 3b), the 
U and I band values deviate significantly from a CCM law, 
which is recovered after the full ew Si and EW Ca correction 
(Eq. lc, Fig. 3c). The empirical fit presented in subsection 3.2 
can be forced to follow a CCM extinction law by substituting 
(SA^, = (cia + bj/Rv) 5A* vj + r\x, where ax and bx are the wave- 
length dependent parameters given in Cardelli et al. (1989) and 
O'Donnell (1994), and a single Ry is fit over all bands. This fit 
applied to SA\ i+Ca leads to an average Ry = 2.78 ± 0.34 for the 
SN host galaxies in our sample. This value is compatible with 
the Milky Way average of Ry = 3.1 The quoted uncertainty is 
statistical and derived with a jackknife procedure, removing one 
supernova at a time. 

We tested the robustness of our Ry determination in several 
ways. Since the sx are measured in sequence, each one using the 
corrected magnitudes from the previous step, we swap the order 
in which the correction is applied, and find Ry = 2.70. Using 
the whole sample to compute the sx instead of applying a 2cr 
clipping cut leads to R v = 2.79, showing the small influence of 
clipping at this stage. Host galaxy subtraction is performed using 
the full spatio-spectral information from the host obtained after 
the SN la has faded, and we do not see evidence for residual host 
galaxy features in our spectra. The measurement covariance ma- 
trix C, is dependent on the assumed calibration accuracy: for the 
present result, this is estimated from repeated measurements of 
standard stars. Another estimate can be obtained from the resid- 
uals to a SALT2 light-curve fit. Using this in C, values leads to 



4 



N. Chotard et al. (Nearby Supernova Factory): Reddening law of SNe la 



Ry = 2.76. The stability of the result with respect to our choice 
of bandpasses is checked by removing one band at a time, and we 
obtain a RMS of +0.23 around the mean result. Finally, a Monte 
Carlo simulation was performed using the converged jx, SA vi , 
r]x and dispersion matrix D as an input, and setting Ry = 2.78. 
The 6A,\j are then randomly generated adding Gaussian noise 
with covariance C,- + D. Over 100 generations, the D matrix is 
recovered with a maximal bias of 5% on the diagonal elements 
and the mean fitted value is Ry = 2.68 ± 0.03, indicating that, if 
anything, our method slightly underestimates Ry. 

5. Discussion 

In this work, we selected ew Si as a first variable for correction 
as it provides a good proxy for X\, and it is a model-independent 
variable. In our data-set the Pearson correlation coefficient of 
ew Si with x\ is -0.83 ± 0.04 (Fig. 1), confirming the result found 
in Arsenijevic et al. (2008). Computing the Hubble residuals in 
the B band using (x\ , c) or (ew Si , c) to standardize SNe both lead 
to a dispersion of residuals of 0.16 mag, also confirming this 
result. We find that ew Si is uncorrected with c(p — -0.04 ±0.1 2, 
Fig. 1), as opposed to the recent claim by Nordin et al. (2010). 
As ew Si and x\ are highly correlated, it is possible to redo the 
analysis with X\ in place of ew Si . The resulting SED correction 
vector s x * is similar to and the conclusions identical, yielding 
R v = 2.69. '' 

Computing the Hubble residuals with (ew Si , EW Ca , c) to stan- 
dardize supernovae leads to a dispersion of 0.15 mag, which is a 
small improvement. Indeed, we do observe a correlation of EW Ca 
and c, with p = 0.34 + 0.10. This correlation is even increased 
to p = 0.50 when c is computed with a U band in addition to 
BVR in the light-curve fit. This shows that c contains an intrin- 
sic component, and that the SALT2 analysis is already account- 
ing for some of this Ca effect. This is illustrated in Fig. 3b: the 
SALT2 color model is derived taking into account a variability 
measured by one intrinsic parameter only, and thus can be com- 
pared to yf. The observed U band contribution in broad filters 
explains the UV rise in the empirical color models. The presence 
of a second intrinsic parameter induces an increased variability 
in the UV as was already noticed by Ellis et al. (2008), but with- 
out a clear attribution to the Ca n H&K line. 

Much lower effective Ry values have been found previously: 
R v = 1.1 in Tripp (1998), R v = 2.2 in Kessler et al. (2009); 
Guy et al. (2010) and R v ~ 1 - 2 in Folatelli et al. (2010). These 
values were derived accounting only for one intrinsic parameter 
beyond color. If we derive an effective value after the sole ew Si 
correction to mimic these analyses, we obtain R v = 3.1, so the 
explanation for the difference has little to do with the number 
of intrinsic parameters entering the correction. This difference is 
explained by the assumption on the dispersion matrix. Indeed, if 
we instead use the one of Guy et al. (2010), corresponding to an 
RMS between 0.09 and 0.11 mag on the diagonal, and all off- 
diagonal terms with an identical RMS of 0.09 mag, we obtain 
Ry = 1.86. This value is representative of previous analyses, but 
the matrix used is incompatible with our best fit matrix. D has 
diagonal values of 0.07, 0.04, 0.05, 0.05 and 0.10 for UBVRI 
respectively, lower than the values quoted by Guy et al. (2010). 
However, while adjacent bands are correlated in our matrix, an 
anti-correlation growing to -1 when the wavelength difference 
increases is observed for other bands, which implies a large color 
dispersion. As shown by Freedman et al. (2009) and Kessler 
et al. (2009), increasing the color dispersion leads to higher R v 
value. This long-range anti-correlation implies that uncorrected 
variability in the SN spectra and/or the reddening law remains. 



6. Conclusions 

Equivalent widths are essentially independent of host reddening 
and provide a handle on the intrinsic properties of the SNe la. In 
particular, ew Si measured at maximum light is confirmed to be 
as powerful as stretch. After correction for this intrinsic contri- 
bution to the magnitude, we have derived a reddening law, and 
proposed a natural method to assess the dispersion due to addi- 
tional fluctuations. We found that this empirical reddening law 
is affected by features linked to the SN physics, and showed that 
the Ca n H&K line provides a second intrinsic variable uncorre- 
cted with ew Si or x\ . Correcting Hubble residuals with ew Si and 
EW Ca leads to a reddening law consistent with a canonical CCM 
extinction law, while adding a dispersion in colors during the fit 
leads to a value of Ry close to the Milky Way value of 3.1. Due 
to the coupling of Ry with the residual dispersion matrix, the de- 
rived value of Ry may be affected as the remaining variability 
becomes better understood and corrected. Our plan to study the 
same SNe la over a range of phases should help address this is- 
sue. The evolution of SNe la will be one of the main limitations 
to their use in precision cosmology surveys at high redshift. Our 
findings show that the accurate measurements of the Ca n H&K 
line provides an additional tool to improve the separation of dust 
and intrinsic variability. The presence of remaining scatter offers 
the possibility of improvements in the future. 
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